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Abstract 

We present a calculation of the n'd scattering length with an accuracy of a few percent using chiral perturbation 
theory. For the first time isospin-violating corrections are included consistently. Using data on pionic deuterium 
and pionic hydrogen atoms, we extract the isoscalar and isovector pion-nucleon scattering lengths and obtain a'^ = 
(7.6 ± 3.1) • 10'^ and a~ = (86.1 ± 0.9) • 10'^ M~^. Via the Goldberger-Miyazawa-Oehme sum rule, this leads to 
a charged-pion-nucleon coupling constant gl/4n: = 13.69 ± 0.20. 

Keywords: Pion-baryon interactions, Chiral Lagrangians, Electromagnetic corrections to strong-interaction 
processes, Mesonic, hyperonic and antiprotonic atoms and molecules. 



1. Introduction 

Hadron-hadron scattering lengths are fundamental 
quantities characterizing the strong interaction, and are 
slowly becoming accessible to ab initio calculations in 
QCD |1, 2]. Among them, of particular interest are 
pion-hadron scattering lengths: the chiral symmetry of 
QCD and the Goldstone-boson nature of the pions dic- 
tate that they are small and their non- vanishing size 
is linked to fundamental quantities like the light quark 
masses and condensates. Chiral symmetry in particu- 
lar predicts that the isoscalar pion-nucleon scattering 
length a'^ is suppressed compared to its isovector coun- 
terpart a~. A precise determination of a'^ would im- 
prove knowledge in many areas, e.g., dispersive anal- 
yses of the pion-nucleon cr-term [4], which measures 
the explicit chiral symmetry breaking in the nucleon 
mass due to up and down quark masses, and is, in turn, 
connected to the strangeness content of the nucleon. 
But, lack of TT^ beams and neutron targets makes di- 
rect pion-nucleon scattering experiments impossible in 
some charge channels, complicating a measurement of 
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the only hope for future access to the Ti^p scatter- 
ing length lies in precision measurements of threshold 
neutral-pion photoproduction Thus, the combina- 
tion of data and theory has, until now, lacked sufficient 
accuracy to even establish definitively that 9^ 0. a", 
on the other hand, serves as a vital input to a determi- 
nation of the pion-nucleon coupling constant via the 
Goldberger-Miyazawa-Oehme (GMO) sum rule 
While the uncertainty in a~ is much less than that in 
a^, it still contributes significantly to the overall error 
bar on the sum-rule evaluation Jvl S] . This is one of 
several examples where data on pion-nucleon scatter- 
ing affects more complicated systems like the nucleon- 
nucleon {NN) interaction, and hence has an impact on 
nuclear physics. 

2. Pionic atoms 

Within the last ten years new information on pion- 
nucleon scattering lengths has become available due to 
high-accuracy measurements of pionic hydrogen (nH). 
The most recent experimental results |9] are 

61, = (-7.120 + 0.012) eV, Yu = (0.823 ± 0.019) eV, 

(1) 

for the (attractive) shift of the 1^ level of nH due to 
strong interactions and its width. These are connected, 

October 29, 2010 



respectively, to the 7r"-proton scattering length, Uj^-p, 
and the charge-exchange scattering length in the same 
channel ifioll . 61^ is related to Uj^-p through an improved 
Deser formula lllll 



(2) 



where a = e^/47i, jjh the reduced mass of tt//, 
K, = 2a(l-\oga)fiHa,-p, and^^ = 2S^h(0)/'¥h(0) = 
0.48% is the effect of vacuum polarization on the wave 
function at the origin B12|] . Further, the width is given 

Tu = 4a'filp,{l + + /^r + C^), (3) 



with 



Kr = 4a(l- \oga)iiHan-p 



(4) 



Here mp, mn, Mj^, and M^jO are the masses of the proton, 
the neutron, and the charged and neutral pions, respec- 
tively, pi is the momentum of the outgoing mfi pair, and 
the Panofsky ratio lUT 



cr{7i p TT^n) 
cr(7i~p ny) 



= 1.546 + 0.009 



(5) 



incorporates the effect due to the radiative decay chan- 
nel of ttH. The pertinent scattering lengths are related 
to a- via (Bl] 



+ a + AcLj^-p, 



«rp = -V2a- + A«-. (6) 



Throughout we follow the notation of Ref. [15] for the 
different nN channels, and have a'^ as a'^ plus a fixed 
shift explained below (see Sect. 13.41) . The other shifts 
in Eq. © take values Aa^^-p = (-2.0 ± 1.3) • 10-^M;^ 
and Aal^J'p = (0.4 ± 0.9) • lO-^M;^ d. This accounts 
for isospin- violating effects up to next-to-leading order 
(NLO) in the chiral expansion. 

Equations 0, ©, and (O permit an extraction of 
a~ and a^. However, further experimental informa- 
tion leads to better control of systematics and could en- 
hance the accuracy of the scattering-length determina- 
tion. Consequently, additional measurements of pion- 
nucleus atoms are of high interest — especially for atoms 
with isoscalar nuclei, as they provide better access to a^. 
Here we use state-of-the-art theory to perform a com- 
bined analysis of the recent data for pionic deuterium 
(nD) as well as the numbers in Eq. ([T]) for nH. The 
resulting values for a~ and a'^ are of unprecedented ac- 
curacy. 
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Figure 1: Topologies for n'd scattering. Solid, dashed, and wiggly 
lines denote nucleons, pions, and photons, respectively. The blobs 
indicate the deuteron wave functions. 



In this work we focus on the strong shift, 6[^, of the 
l^" level of pionic deuterium, which is related to the real 
part of the 7r"-deuteron scattering length. Re aj^-d, by an 
improved Deser formula analogous to Eq. ^ [ I16 [| : 



-2a'filRca,-d{l^KD^Sll'). (7) 



In Eq. O we have = 2^^z)(0)/^d(0) = 0.51% 
I I2I] , Kd = 2a(l - loga)yUz)Rea;r-j, and jid as the nD 
reduced mass. 

3. The pion-deuteron scattering length 

The real part of Gj^-d can be decomposed into its two- 
and three-body contributions as: 



Rca,-d = af}^^a% 



(8) 



It is in a^^}r that resides. Therefore, a^^^r must be 

7T d ^ 7T d 

calculated reliably if measurements of are going to 
be profitably exploited to get information on . 

Thus, the bulk of the rest of this paper describes a 
calculation of chiral perturbation theory OifPT). 

This quantity can be expressed as 



(3) 

7T d 



(9) 



where a^^^ defines the strong contribution, ^^^^p^^ in- 
volves two-nucleon and A-isobar-nucleon intermediate 
states, as well as diagrams with crossed pion lines, and 
a™ involves photon-exchange contributions. This last 
piece is present because isospin violation from the up- 
down quark mass difference and electromagnetic effects 
must be taken into account (as in Ref. fl^ we use a 
counting where e - p). Consistent consideration of 
such effects is a key advance made in this paper. We 
now deal with each of the contributions in Eq. ©, be- 
fore returning to a^^? ,. 

3.1. Strong contributions (a^^^) 

The leading diagrams contributing to a^^^ are shown 
in the first line of Fig. [T] So far no counting scheme 
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is known that permits consistent, realistic, and simulta- 
neous consideration of the two- and three-body opera- 
tors which contribute to 7i~d scattering. However, each 
of these operators can be calculated independently, i.e. 
within its class, with a controlled uncertainty. In par- 
ticular, Ref . IIITIi showed how the original counting by 
Weinberg can be modified such that the three- 

body contributions to Uj^-d are calculated to very high 
accuracy. Since isospin breaking in the two-body sec- 
tor is also well under control llisl] . this permits a precise 
extraction of . Therefore, we now discuss the power 
counting for all contributions to a^^^ relative to the lead- 
ing, (9(1), diagram (Ji). 

In this counting there is a (N^N)^7i^7i contact term as- 
sociated with the short-distance pieces of the integrals, 
which enters with an unknown coefficient at O(p^). This 
contribution cannot easily be determined from data, and 
is a key source of uncertainty in our result. With p ~ 
Mj^/ntp, we anticipate an accuracy of a few per cent for 
threshold n'd scattering. This expectation is substanti- 
ated by the sensitivity of our integrals to the choice of 
the deuteron wave function (see below). There we see a 
residual scale dependence of about 5%: an independent 
estimate of the contact term's effect. 

But, to reach this accuracy, we must include all three- 
body terms up to 0(p^^^). In Ref. [59] it was shown that 
the sum of all NLO, 0(p), contributions vanishes in the 
isospin limit, corrections to which only enter at O(p^). 
Thus, the diagrams we need to consider up to 0(p) are 
(di)-(d4) in Fig. [T] Note that although we count (J5) 
as O(p^), its value is enhanced by a factor of tt^ due 
to its topology of two successive Coulombic propaga- 
tors llivl boh . Similar enhancements are present for all 
terms of the multiple- scattering series. Despite this, the 
multiple- scattering series converges quite quickly: we 
find from an explicit calculation that the sum of the first 
two terms (di) and (J5) differs from the full result by 
only 0.1 • 10~^M~^. Note that the next diagram, where 
the pion leaves the two-nucleon system after four nN 
interactions on alternating nucleons, is logarithmically 
divergent, and therefore seems to necessitate a contact 
term. As the terms in the multiple- scattering series are 
enhanced as just described, we expect this contact term 
to also be enhanced. However, that enhancement is not 
enough to overcome the p^ suppression relative to the 
leading, double- scattering, piece of Gj^-d, and so any 
such contact term has an appreciably smaller effect than 
the O(p^) contact term. Therefore its contribution does 
not impact the uncertainty estimate given above. 

To achieve the requisite accuracy for our extrac- 
tion we also need to include isospin-violating correc- 
tions from the different masses of the proton and neu- 



tron and charged and neutral pions in the diagrams (di)- 
(J4). We then express the sum of diagrams in the first 
row of Fig. [T] as: 



(10) 



The first four terms arise from diagrams (di) and (J2). 
However, (J2) is partly accounted for in the two-body 
contribution a^^]^. In order to treat the three-body dy- 
namics properly we must replace the contribution of 
the two-bo dy ( nN) cut there by that of the three-body 
(ttNN) cut MM- The necessary integrals can be rear- 
ranged as in Eq. ([TOl) (for details see 11221] ). a^^^^^^ cor- 
responds to (di) evaluated with a static pion propaga- 
tor, and is numerically by far the dominant contribution. 
a^Lo incorporates recoil corrections to the static pion 
propagator; a^^^ comprises effects due to the three-body 
n^nn and n'pn cuts, and Aa^^^ emerges as an isospin- 
violating correction in this rearrangement. (In principle, 
there are also contributions with P-wave interactions be- 
tween nucleons in the intermediate state, but they are of 
higher order.) Finally, a^^ in Eq. ([TOl) is determined by 
(J3) and (J4), while a^^^P^^ results from (^5). Isospin- 
breaking corrections to the nN scattering lengths that 
appear in a^^^ are relevant only for a^^^^^^, to which they 
contribute about 1%. 

Our power counting is based on dimensional analysis 
assuming all integrals scale only with Mj^. In fact, the 
integrals in Eq. ([TOl) involve other scales too: ^[M^ — 
due to the three-body cut — and ^^/m^, thanks to the 
deuterium wave functions (e is the deuteron binding en- 
ergy). At first glance, the presence of a three-body cut in 
the integral for a^^^ makes it appear to be enhanced over 



its naive ;^PT order by ^Jm^JW^ ll23[l . However, this 
turns out not to be the case, because the Pauli principle 
and the spin-isospin character of the leading nN scat- 
tering operator ensure that the intermediate NN state in 
{di) -h (J2) is projected onto a P-wave 112 111 . In conse- 



do not enter the fi- 



quence the scales ^Mj^e and 
nal result: any enhanced contribution cancels due to a 
subtle interplay between the two diagrams that is dic- 
tated by the Pauli principle. The combined integral is, 
as originally assumed in establishing the ;^PT ordering 
of diagrams, then dominated by momenta of order Mj^. 

The results for the pieces of a^^^ are given in Table [T] 
They produce a total: 



= (-22.6+1.1+0.4) 



10-^m:^ 



(11) 



The first error comes from the evaluation of all men- 
tioned diagrams using different deuteron wave functions 
(we use NNLO chiral (five wave functions with differ- 
ent cutoffs) H, CD Bonn Q, and AV 18 S poten- 
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Table 1: Strong contributions to a^^l. in units of 10 ^M^^ . Here and 
below results are quoted for a~ - 86.1 • 10~^M~^. For the band in 
Fig.[2]the full a' dependence is taken into account. 



^static 


-24.1 + 0.7 


^static 
^NLO 


3.8 + 0.2 




-4.8 + 0.5 


^triple 


2.6 + 0.5 


(3^^ 


-0.2 + 0.3 


Aa(2) 


0.2 



tials), while the second is due to the uncertainty in the 
isospin-breaking shifts in the nN scattering lengths lfl5l] . 

3.2. Photon loops (a^^) 

Effects in this class due to photons with momenta of 
order aMjj; are included in observables via the improved 
Deser formula. Thus, our calculation of a^^}, should 
include contributions from momenta above aMj^. The 
leading contributions due to the exchange of (Coulomb) 
photons of momenta of order Mj^ between the tt" and the 
proton are shown in the second row of Fig. [T] {d^), (dj), 
and (Jg). Photon exchange is perturbative at |k| ~ Mj^ 
(in contrast to the hadronic-atom regime where the pho- 
ton ladder needs to be resummed), and the pertinent 
pieces of these graphs enter at 0(p) relative to (di). 
Such effects in the other diagrams are of a higher ;^PT 
order than we are considering here. 

However, diagrams (de) and (d^)-(dio) are reducible 
in the sense originally defined by Weinberg lll8ll . with 
the ttNN intermediate state involving relative momenta 
of order sjMj^e <^ Mj^. Furthermore, in these diagrams, 
this state can occur with the NN pair in an S -wave, so 
we must also allow for the possibility of NN interactions 
while the pion is "in flight". When this is done we see 
that these four diagrams have an infrared divergence in 
the limit 6^0, being enhanced by ^Mj^je as compared 
to their naive ;^PT order. 

In order to avoid double counting we must also 
subtract from the resulting expressions for {d^ and 
{d%)-{d\{^ (plus NN intermediate- state interactions) 
the quantum-mechanical interference between a zero- 
range (strong) pion-deuteron potential, proportional to 
Qj^-d, and the Coulomb interaction. That interference 
is already accounted for in the improved Deser for- 
mula d?]). Note though, that Eq. d?]) only accounts for 
intermediate- state pion (and deuteron) momenta of or- 
der aMj^. In particular, deuteron structure plays no role 
in its derivation. 

After the pieces of {d^ and {d%)-{d\^ that are already 
included in Eq. d?]) are removed the result is finite. The 
remaining, finite parts of {d^ and {d%)-{d\{^ capture the 
effects of momenta » aMj^- in these loops. These con- 
tributions are defined here to be part of a^-d, and must 



be calculated explicitly. In particular, they include ef- 
fects in the loop which arise from the electromagnetic 
and pion-deuteron "form factors": the manner in which 
the finite extent of the deuteron modifies the loop inte- 
gral for momenta well above the hadronic atom scale 

This contribution to aj^-d is ostensibly large, since it 
is an infrared- sensitive integral that potentially has con- 
tributions from momenta of order ^Mj^e. But, analysis 
analogous to Ref. |0] shows that this particular piece 
of the integral is zero because of symmetry arguments. 
When the NN pair is in an S -wave it can be written as a 
sum of overlaps between NN wave functions in the con- 
tinuum and the deuteron bound state, and orthogonality 
then guarantees that the result is zero. In the case of 
an intermediate NN P-wave pair it is the Pauli principle 
that causes the cancellation [22] . 

There is still a possible contribution in the loop from 
momenta of order ym^. This would be enhanced by 
^nl V^^^ compared to its naive ;^PT order, and so 
could be relevant for our analysis. Direct evaluation of 
this part of {d^ and (Jg) — including the diagrams with 
A^A^ interactions in an S -wave — yields a contribution to 
ajj^-d of -0.04 a^. (Isospin-breaking shifts of a" can be 
added here, but do not change the pref actor.) Replacing 
the single nN scattering of these diagrams by double 
scattering as in (J9) and (Jio) gives effects larger by a 
factor of a^-dl^a" ■> but, despite their being infrared en- 
hanced, the impact of such pieces on aj^^-d is still signif- 
icantly less than our theoretical uncertainty. 

This leaves us needing to consider only effects from 
momenta M;^ in diagrams {d^-{di). As with {di) in 
parts of these diagrams are already included in cFl^, but 
this can be dealt with along the same lines \Tm . The 
result is: 



(P^ = (0.94 + 0.01) • io-'m;\ 



(12) 



where the error again reflects the wave-function depen- 
dence. Thus, virtual photons with |k| ~ M;^ increase 
Re Gj^-d by about 4%. 

3.3. Dispersive and Delta( 1232 ) corrections ( a^^^P^^ j 

These produce effects in a^^-d that scale with half- 
integer powers of p [28, 29]. Their leading contribution 
is 0(p^^^) relative to (Ji), and is computed here using 
a calculation for NN dn up to NLO in ;^PT f^. 
Note that although we include Delta(1232) effects in the 
ttA^A^ TiNN transition operator, it is not necessary to 
account for the Delta(1232) as an explicit degree of free- 
dom when computing the deuteron wave function. Its 
effects in the A^A^ potential at energies of order e enter 
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Figure 2: Combined constraints in the cL"-a plane from data on the 
width and energy shift of nH, as well as the nD energy shift. 



only at relative 0{p^) iQ. In Refs. all integrals 

were cut off at 1 GeV; we have checked that this does 
not introduce additional uncertainty and obtain: 



a^^^P^^ = (-0.6±1.5)- 



10-^m:^ 



(13) 



Since this is at the limit of our desired accuracy we need 
not include isospin- violating corrections to a^^^P^^. 

3.4. The two-body part (ci\-J 

As alluded to above, it is not possible to isolate a" 
in analyses of nH and nD. Information on the isoscalar 
scattering length can only be extracted as a combination 
a^, in which the low-energy constants ci (which occurs 
because its impact on a" is proportional to the neutral- 
pion mass squared) and fi (which denotes the leading 
isoscalar electromagnetic correction) also appear 



1 



Ml 



ci - 2a fi \. (14) 



1 + MJm^ I TiFl 



In the two-body part of Uj^-d, is further shifted, as 
shown in the NLO analysis of Ref . ifish : 



(2) 
n a 



{a + Aa ) , 



Aa-" = (-3.3 ± 0.3) • lO'^M;^ 



(15) 



4. Results and Discussion 

We now add together all the individual contribu- 
tions. Amusingly, most of the additional three-body 
corrections considered in this study accidentally cancel: 
Aa(2) + ^static _^ ^cut _^ ^EM = (0.1 + 0.7) • 10-^M;1. For 



Table 2: Individual contributions to the error on are added in 
quadrature to obtain the uncertainty depicted in the bands of Fig. |2] 
Each row below gives the impact of one source of error as a percent- 
age of that total. The first row is the impact of the experimental uncer- 
tainty in e[^, the second gives the uncertainty in the isospin-breaking 
shifts of nN scattering lengths that occur in a^^^, and the third row is 
the uncertainty in Aa'^ according to Eq. ( 115) . The final two rows show 
the impact of uncertainties in our calculation of Rea^^J^, as described 
in the text. 





16% 


Aa^-p.Ad"^^^ 


21 % 


Act 


30% 


^disp+A 


75% 


Wave-function averages 


53% 



this reason, the main impact of our analysis on the ex- 
traction of pion-nucleon scattering lengths turns out to 
be due to the NLO isospin-breaking corrections in the 
two-body part ll3lll . 

The energy shift of nD has recently been remeasured 

as m 

6f^ = (2.356 + 0.031) eV. (16) 

Combining this result, the dependence of the n~d scat- 
tering length on a" and a~ , and the results for nH dis- 
cussed above, we find the constraints depicted in Fig. [J] 
The combined Icr error ellipse yields 

a - = (1.9 + 0.8)- 10"^M;^ a- = (86.1 +0.9)- 10"^M;^ 

(17) 

with a correlation coefficient pa-a+ = -0.21. We find 
that the inclusion of the nD energy shift reduces the un- 
certainty of a'^ by more than a factor of 2. Note that 
in the case of the nH level shift the width of the band 
is dominated by the theoretical uncertainty in Aa^-p, 
whereas for the nH width the experimental error is 
about 50% larger than the theoretical one. The uncer- 
tainty in a^^^P^^ is the largest contribution to the nD er- 
ror band, see Table [2l The wave-function averages con- 
tribute about 0.5 • 10"^M^^ to the overall uncertainty in 
a" , which is in line with the estimated impact on a^^-^ of 
the 0{p^) — relative to {d\) — contact term. 

Taken together with ci = (-1.0 + 03) GeV" ^ Q and 
the rough estimate |/i | < 1 .4 GeV"^ I!! , Eq. ([17]) yields 
a non-zero a" at better than the 95% confidence level: 



a - = (7.6 + 3.1) 



10-^m:^ 



(18) 



The final result for is only a little larger than 
several of the contributions considered in our analysis. 
This emphasizes the importance of a systematic order- 
ing scheme, and a careful treatment of isospin violation 
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and three-body dynamics. A reduction of the theoreti- 
cal uncertainty beyond that of the present analysis will 
be hard to achieve without additional QCD input that 
helps pin down the unknown contact-term contributions 
in both the nN and nNN sectors. 

Finally, these results allow us to infer the charged- 
pion-nucleon coupling constant, gc, from the GMO sum 
rule, with isospin- violating corrections to the nN scat- 
tering lengths fully under control for the first time. In- 
spired by Ref. JtI], we take Uj^-p extracted from Eq. ©, 
aj:-p + Gj^-n from our Gj^-d analysis, and Gj^-n - Gj^+p from 
Ref. Q, yielding g^/An = 13.69 ± 0.12 ± 0.15 iQ. 
(Here the first error is due to the scattering lengths and 
the second to an integral over n-p cross sections jvilsll.) 
This is i n ag reement with determinations from NN ll34ll 
and TiN 11351] scattering data. 
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